A Layer-Stripping Solution of the Inverse Problem for a One-Dimensional Elastic Medium

A fast algorithm for recovering profiles of density and Lame parameters as functions of depth for the inverse seismic problem in an elastic medium is obtained. The medium is probed with planar impulsive P and SV waves at oblique incidence, and the medium velocity components are measured at the surface. The interconversion of P and SV waves defines reflection coefficients from which the medium parameter profiles are obtained recursively. The algorithm works on a layer-stripping principle, and is specified in both differential and recursive forms. A physical interpretation of this procedure is given in terms of a lattice filter, where the first reflections of the doangoing waves in each layer yield the various reflection coefficients for that layer. A computer run of the algorithm on the synthetic impulsive plane wave responses of a twentylayer medium shows that the algorithm -woriks satisfactorily. *The work of this author was supported by the Exxon Education Foundation. **The work of this author was supported by the Air Force Office of Scientific Research under Grant AFOSR-82-0135A.